Spin diffusion of correlated two-spin states in a dielectric crystal 
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Reciprocal space measurements of spin diffusion in a single crystal of calcium fluoride (CaF2) 
have been extended to dipolar ordered states. The experimental results for the component of the 
spin diffusion parallel to the external field are = 29±3x 10"^^ cm /s for the [001] direction and 
Di!) = 33±4x 10-^2 cm /s for the [111] direction. The measured diffusion rates for dipolar order are 
faster than those for Zeeman order and are considerably faster than predicted by simple theoretical 
models. It is suggested that constructive interference in the transport of the two spin state is 
responsible for this enhancement. As expected the anisotropy in the diffusion rates is observed to 
be significantly less for dipolar order compared to the Zeeman case. 

PACS numbers: 75.45.+j,76.60.-k 



In this work we measure and compare the spin diffu- 
sion rates for dipolar and Zeeman energy in a dielectric 
crystal. Spin diffusion provides a well-posed problem in 
the study of multi-body dynamics, as the Hamiltonian 
of the system is well known. The study of the diffusion 
of the Zeeman and dipolar ordered states is essentially 
the study of the evolution of different initial states un- 
der the secular dipolar Hamiltonian. In particular the 
Zeeman ordered state consists of single spin population 
terms only, while the dipolar ordered state consists of 
correlated two spin states In a strong magnetic field 
these quantities are independently conserved, and have 
different diffusion rates and spin-lattice relaxation times. 
The spins are in dynamical equilibrium under the action 
of the secular dipolar Hamiltonian, and the constants of 
the motion are only defined for the total spin system. 
If the identity of individual spins can be distinguished, 
the spins are clearly seen to evolve in time as in the case 
where we write a position dependent phase on the spins. 
This is the essence of our spin diffusion experiments. Pre- 
vious attempts at measuring spin diffusion of Zeeman 
0, 1^ and dipolar energy relied on theoretical mod- 
els that related the diffusion rate to observed relaxation 
times. Zhang and Cory used reciprocal space techniques 
to perform the first direct experimental measurement of 
spin diffusion in a crystal of CaF2 0. 

There have been a number of attempts to calculate the 
rate of spin diffusion in single crystals, for both the Zee- 
man energy and the dipolar ener^ of the spin system 
using theoretical models 0, S |l E [H [H [H , and 
classical simulation 0, 0| . These calculations suggest 
that the spin diffusion rates of Zeeman and dipolar order 
should approximately be the same. However, an exami- 
nation of the physical processes leading to the diffusion 
shows that the dynamics can be quite different. Figure 1 
shows a simple illustrative model that suggests that the 
diffusion of dipolar order should be faster (and more com- 
plicated) than that of Zeeman order, due to the increase 
in the number of possible paths for the propogation of 
the dipolar ordered state. While a faster diffusion rate 
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FIG. 1: (a) For the diffusion of Zeeman order along the z 
direction spins 1 and 2 need to undergo a flip-flop, (b) For 
the diffusion of dipolar order along the z direction, there are 
two possibilities, both spins 1 and 3, and spins 2 and 4 can flip, 
or spins 1 and 2, and spins 3 and 4 can flip. There are thus 
two different paths to the same final state and interference 
effects may be observed, (c) If spins 1 and 2 are initially in 
the same state (both | or both j), no evolution takes place, 
(d) Even if states 3 and 4 are initially ] and [, two different 
evolution paths are present. Spins 1 and 2 can flip and spins 
3 and 4 can flip or spins 1 and 4 and spins 2 and 3 can flip. 
Thus dipolar diffusion dynamics are less easily quenched. 



would seem to make the measurement of dipolar diffu- 
sion easier, the measurement is more challenging as the 
experiment has to be carried out in a time comparable 
to Tio (usually an order of magnitude smaller than Ti). 

The basic reciprocal space NMR spin diffusion experi- 
ments for Zeeman and dipolar order are outlined below. 
The Hamiltonian of the system is given hyU. — Hz+Hd, 
where Tiz = huj ^ ^ is the Zeeman interaction, and 
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is the dipolar interaction, lu is the Larmor frequency of 
the spins and 6^ = 7_ji i-3cos e,, ^j^^ strength of the 

dipolar coupling between spins i and j. In the high field 
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(|Wz| >> I'WdI), high temperature {f3huj << 1) limit, 
the initial equilibrium density matrix is given by 



Po 



1 (^l~f3f^Y.r}j=^-6p 
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where Z is the partition function and 1 is the identity. 
In the following discussion we follow the evolution of 6 p. 

After the (7r/2)y pulse, the magnetization grating is 
written onto the spin system by interleaving strong mag- 
netic field gradient pulses in the long delay times within 
a magic echo train 0, 0| . The magic echo train sus- 
pends evolution of the dipolar coupling, and the effective 
Hamiltonian of the system is just the gradient term. 



(3) 



where 7 is the nuclear gyromagnetic ratio, dB^/dz is the 
gradient field strength, and Zi is the spatial location of 
spin i. This creates spatially modulated Zeeman mag- 
netization. If the total time the gradient is applied is 
(5, the wave number corresponding to the applied gradi- 
ent is fc = 7(5 [dBz/dz). Following the application of the 
magnetic field gradient, the state is 



5p 
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where (f>i = kzi tags the spatial location of spin i. The 
state in Eq. (4) is no longer a constant of the motion 
under a collective measurement of the spins. A ['K/2)y 
pulse flips one component of the transverse magnetiza- 
tion back along the longitudinal (2) axis. The remaining 
transverse component decays rapidly yielding 
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These steps are common to both experiments. 
Zeeman diffusion 

We allow the above state to evolve under Hd for a time 
A. Following spin diffusion, we can express the state as 
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where the coefficients Cij(A) represent the relative 
amount of polarization initially at site i that has been 
transported to site j after time A. We have also ne- 
glected higher order spin terms in the above equation, as 
these do not contribute when we trace over the spins. In 
the absence of the phase tag, these terms would sum to 
zero, as the Zeeman state is a constant of the motion. 

In the final step, a (7r/2)y pulse is applied followed by a 
magnetic field gradient, equal in amplitude but reversed 
in direction to the encoding step. This gradient labels 



the new location of the polarization, in analogy with a 
scattering experiment. Following a (7r/2)y pulse, the lon- 
gitudinal component of the magnetization is 



SpiA) = ^ 



mi 
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In the absence of spin diffusion, Cy = (5^ , the phase terms 
would be identical and the original state would be recov- 
ered (scaled by 1/2). In the presence of spin diffusion 
the residual phase cos((/)i — (pj) reflects the spatial trans- 
port of the magnetization. In the experiment, the total 
z magnetization is measured, tracing over the spatial co- 
ordinates of the sample. Under diffusive dynamics the 
observed attenuation is exp(— A;^Z3A). 
Dipolar diffusion experiment 

The dipolar diffusion experiment requires two addi- 
tional steps. The modulated Zeeman magnetization (Eq. 
(5)) is first converted to encoded dipolar order by ap- 
plying either the two-pulse Jeener-Broekaert (JB) se- 
quence |l8| | or an adiabatic demagnetization in the rotat- 
ing frame ( ADRF) . Following the creation of dipolar 
order, the state of the system is 

This evolves under T-Ld for time A during which spin 
diffusion takes place. The state can then be expressed as 



5p 



I3h 



i<j;k<l 



(A) cos 0i X 



-{iliL+iiii: 



(9) 



We have again left out higher order terms that do not 
contribute to the measured magnetization. These terms 
would sum to zero in the absence of a phase tag, as 
the dipolar ordered state is a constant of the motion. 
The two-spin dipolar ordered state is converted back to 
a single-spin Zeeman state using either a 7r/4 pulse or an 
adiabatic remagnetization in the rotating frame (ARRF) . 
The 7r/4 pulse only converts a portion of the dipolar or- 
dered state back to Zeeman order, while the ARRF allows 
a complete refocusing to the Zeeman state. The gradient 
modulation is subsequently unwound yielding 



6p- 
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(10) 

where eki is a dimensionless quantity proportional to bki ■ 
In the limit of infinitesimal polarization the dipolar or- 
dered states created are local two-spin correlations. The 
length scale of these correlations (A) is much smaller than 
the length scale of the spatial modulation (lOO's of nm), 
yielding 0^ = . As a consequence the phase spread due 



3 




(b) 

0' 

-0.5 
-1 

-1.5 
-2 



k^A (s/cm^^ 



0.5 



2 ■'■5 2 
\C\ (s/cmT 



2.5 



X 10 





□ k/27t=353 1/mm 
^ k/27t=644 1/mm 
> k/27t=1053 1/mm 
o k/27t=1500 1/mm 


1 









(a) 



h-0.5 

[ -1 

-1.5 

( 

(b) 









> 


k/2jt = 636 1/mn 


[ADRF] 








□ 


k/2ji = 659 1/mr 


n [ADRF] 








o 


k/2jt = 943 1/mn 


[JB] 




^ k f r~ 


■ 




k/27r = 991 1/mn 


[JB] 















0.5 



1.5 2 

k^A (s/cm^) 



2.5 



3.5 



xlO" 



^ k/2!i = 461 1/mm[JB] 

O k/2!i = 632 1/mm [ADRF] 

<] k/2ji = 659 1/mm [JB] 

□ k/2ji= 1055 1/mm ]JB] 



X 10 




1.5 2 

k^A (s/cm^ 



FIG. 2: (a) Experimental results of the spin diffusion rate of 
Zeeman order with the crystal aligned along the [111] direc- 
tion. The pulsed gradient fields generated are o 2886 G/cm, 
* 5661 G/cm, < 8510 G/cm, □ 8510 G/cm. (b) Experimental 
results of the spin diffusion rate of Zeeman order with the 
crystal aligned along the [001] direction. The pulsed gradi- 
ent fields generated are □ 2886 G/cm, * 5661 G/cm, > 8695 
G/cm, o 7474 G/cm. The error in k is approximately 3 %. 



FIG. 3: (a) Experimental results of the spin diffusion rate of 
dipolar order with the crystal aligned along the [001] direc- 
tion. The pulsed gradient fields generated are > 5698 G/cm, 

□ 5698 G/cm, o 8473 G/cm , * 8399 G/cm. (b) Experimen- 
tal results of the spin diffusion rate of dipolar order with the 
crystal aligned along the [111] direction. The pulsed gradient 
fields generated are * 2886 G/cm, o 8436 G/cm, < 5661 G/cm, 

□ 8362 G/cm. The error in k is approximately 3 %. 



to the creation and refocusing of dipolar order can also 
be safely ignored in the experiments here. In the absence 
of spin diffusion, (pi = (f>k and the original state is recov- 
ered (scaled by 1/2). As in the Zeeman case, the residual 
phase term cos((/)i — (j)k) encodes the spatial transport in 
the presence of spin diffusion. 

We performed the measurements at room temperature 
at 4.7 T (i^F 188.35 MHz) using a Bruker Avance spec- 
trometer, and a home-built z-gradient probe. The gra- 
dient coil had an inductance of 26 fiK, a resistance of 
1.6 n, and an efficiency of 370 ± 10.5 G/cmA The 
solenoidal RF coil geometry used has been shown to pro- 
duce a homogeneous field over a small sample volume 
|2l| . The experiments were performed on a 1 mm^ single 
crystal of CaPa having Ti = 256.2 ± 3.4 s along the [001] 
direction and Ti ~ 287.5 ± 7 s along the [111] direction. 
The measured Tjo of our crystal was 9.4 ± 0.4 s along 
the [001] direction and 9.3 ± O.I s along the [111] direc- 
tion. The 7r/2 pulse was 1.65 /iS, the spacing between 
magic echoes was 60 /iS and the length of each gradi- 
ent pulse was 35 ^s. The coil constant was determined 
by calibration with a sample of water whose molecular 
diffusion coefhcient is well known. A phase alternating 
scheme was applied to correct for RF pulse errors in the 
magic-echo sequences ,17 j ■ A 7r/4 pulsed refocusing of 
the dipolar ordered state was used in the JB experiment 
and an ARRF was used with the ADRF. 

The measured signal at long times, A ^ (271/6™"^), is 

S ^ Soeiq){-k^DA)eM-NTix) (11) 



TABLE I: Summary of the experimental results of the spin 
diffusion rate of spin-spin energy, Dd, and Zeeman energy, 
Dz for single crystal calcium fluoride. 



This measurement 


[001] 


[111] 


Di^Qx 1 D\\\ 


D^l, (xlQ-^^cmVs) 


29 ± 3 


33 ± 4 


0.88 ± 0.14 


(xlO-i^cmVs) 


6.4 ± 0.9 


4.4 ± 0.5 


1.45 ± 0.26 


Prior measurement \5] 


[001] 


[111] 


Diiw 1 D\\\ 


d\ (xlO'^^cmVs) 


7.1 ± 0.5 


5.3 ± 0.3 


1.34 ± 0.12 


Theoretical studies of d\ 


[001] 


[111] 


DciQx 1 D\\\ 


Ref. [11] (xlO"^^cmVs) 


6.98 


4.98 


1.4 


Ref. [12j (xlO"^^cmVs) 


8.22 


6.71 


1.22 


Ref. yj (xlQ-^^cmVs) 


7.43 






Theoretical studies of D^i, 


[001] 


[111] 


DciQx 1 D\\\ 


Ref- ^ (xlO"^^cmVs) 


8.53 


8.73 


0.98 


Ref. fl4] (xlO"^^cmVs) 


13.3 






Ratio of Dd to Dz 


[001] 


[111] 




This measurement 


4.5 ±0.8 


7.5 ±1.3 




Ref fl3] 


1.22 


1.75 




Ref. [14] 


1.79 







as the system also undergoes spin-lattice relaxation dur- 
ing the time A (^\x — T\ in the Zeeman experiment and 
T\r^ — T\D va. the dipolar experiment). It is necessary to 
correct the signal amplitude for this attenuation in order 
to obtain the correct diffusion coefhcient. The crystal 
used in this experiment has a much longer T\ than the 
crystal used previously for the Zeeman diffusion mea- 
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surement (Ti = 114.2 ± 5.3 s along the [001] direction 
and Ti = 156.8 ± 9.7 s along the [111] direction), since 
the dipolar diffusion experiment needs to be carried out 
in a time comparable to Tio {Ti£, < Ti). In figures 2 
and 3 we plot In(S') + A/Tia; versus fc^A for the Zeeman 
and dipolar experiments when the external magnetic field 
is applied along the [001] and [111] crystal orientations. 
The diffusive character of the dynamics is shown by the 
universal scaling observed with the different gratings. 

Table 1 summarizes the experimental results and shows 
previous theoretical and numerical results for compari- 
son. The measurements and theory agree reasonably well 
for Zeeman order, both in the absolute value and the ori- 
entational dependence of the spin diffusion rate. Though 
the various theories predict the ratio of spin diffusion 
rates for dipolar and Zeeman order should be between 
one and two, the measurement yields a ratio larger than 
four. In contrast to Zeeman diffusion, dipolar diffusion 
exhibits a much smaller anisotropy. This is expected as 
the anisotropy of the field generated by a single magnetic 
dipolc is larger than that created by a pair of dipoles. 

In the description of the Zeeman and dipolar diffusion 
experiments above, the details of the dynamics have been 
aggregated and included in the coefficients Cy and dyfei- 
Evaluation of these terms would entail solving the many 
body problem for spin diffusion. However, in order to 
compare the two diffusion processes, it is illustrative to 
evaluate these coefficients in the short time limit. Eval- 
uating the second commutator for the persistence of the 
initial states, we obtain 

c»'^W = i-7E(^^^)' + o(^') (12) 

d.,,,{t) « 1 _ ^ ^ (62^ + h]^ + huh.u) + 0{t^). (13) 

The equation for dijij above is not exact, as there are 
also some two spin terms present that are not propor- 
tional to the dipolar ordered state. If the dipolar ordered 
state had been proportional to Illl, then we would have 
diju = CikCji- The coefficient d in Eq. (13) differs from 
a simple product of the c's at this order by the cross- 
terms bikbjk- These cross terms could either increase or 
decrease the decay rate as the dipolar couplings can be 
both positive and negative. Thus both constructive and 
destructive interference effects could be obtained in prin- 
ciple. This suggests that the rapid diffusion rate observed 
for dipolar order is a consequence of a constructive inter- 
ference effect that enhances the diffusion rate over that of 
unentangled spin pairs. Currently available calculations 
do not take this interference effect into account properly. 
For instance, the moment method |l2l | introduces an ad- 
hoc cutoff function in frequency space to force an expo- 
nential solution, even though the coherent signal evolu- 
tion should actually be an even function of time. Indeed, 



the cross terms in Eq. (13) cancel out the level of the sec- 
ond moment, and it is unclear whether they remain in the 
higher moments, as the calculation of higher order mo- 
ments has historically proved challenging. Incorporating 
such interference effects in the calculation is important 
in our system where the dynamics remain coherent and 
reversible over the timescale of the experiment - about 
10^ times the correlation time of the dipolar coulping. 

In conclusion, we have performed a direct measure- 
ment of spin diffusion of dipolar energy in a single crys- 
tal of CaF2. The rate is observed to be significantly 
faster than expected from previous theoretical studies, 
though the reduced anisotropy is well predicted. Given 
the highly mixed nature of the initial density matrix in 
these experiments and the large Hilbert space, it is diffi- 
cult to characterize the dynamics in more detail. Measur- 
ing spin diffusion as a function of increasing polarization 
would allow us to follow the spin dynamics as they are 
constrained to progressively smaller subspaces, enabling 
a more careful exploration of the Hilbert space dynamics. 
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